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Black hole perturbation theory, or more generally, perturbation theory on a 
Schwarzschild bockground, has been applied in several contexts, but usually un- 
der the simplifying assumption that the ADM momentum vanishes, namely, that 
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the evolution is carried out and observed in the "center of momentum frame". In 



this paper we consider some consequences of the inclusion of a non vanishing ADM 
momentum in the initial data. We first provide a justification for the validity of 
the transformation of the initial data to the "center of momentum frame" , and then 
analyze the effect of this transformation on the gravitational wave amplitude. The 
most significant result is the possibility of a type of gravitational memory effect that 
appears to have no simple relation with the well known Christodoulou effect. 

> 

^ ! PACS numbers: 04.30.-w 

00 

o 

\o 
o 

O I- INTRODUCTION 

Following its successful application to the close limit of head on collisions, by Price and 
Pullin [1] , black hole perturbation theory has become an important tool in the analysis of the 
final stages of the coalescence of two black holes, after their collision or final plunge from their 
innermost stable orbit Q. In this context, the theory, originally formulated in the Fourier 
^ ■ transformed frequency domain by Regge and Wheeler [3] and Zerilli pj, is considered in the 
time domain as a manner of (perturbatively) solving an initial value problem. Namely, one 
is given initial data (a solution of the constraint equations of General Relativity) in the form 
of the 3-metric and extrinsic curvature on (some region of) a three dimensional spacelike 
hypersurface £ and the problem is to find the full metric in the domain of dependence of the 
initial data. The type of data considered here depends in general on one or more perturbation 
parameters in such a way that one recovers that for a Schwarzschild black hole when the e; 
vanish. The possibility of a perturbative analysis of the evolution is based on the expansion 
of the data in an appropriately chosen angular basis (tensor spherical harmonics), and the 
assumption that the angular component coefficient functions can be further expanded in 
powers of the e,. The non trivial part of the Einstein equations can then be cast in the 
form of an infinite set of coupled partial differential equations for functions of two variables 
(conventionally, the Schwarzschild coordinates t and r), that can, in principle, be integrated 
order by order in the perturbation parameters 0. This, however, does not take into account 
the invariance of the geometry under coordinate changes. In fact, when this invariance is 
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fully accounted for, the relevant physical information ends up being encoded in two sets of 
functions, where each element corresponds to an angular mode, introduced respectively by 
Regge and Wheeler and by Zerilli 0], that satisfy wavelike equations in t, r. 

While the above analysis, leading to the emergence of the Regge - Wheeler and the Zerilli 
functions is based on the general invariance of the theory under coordinate transformations, 
in recent applications to the black hole binary problem one is interested in the evolution of 
a given initial data, posed on a given hypersurface S. In perturbation theory, a particular 
choice is made of the zeroth order (in the e«) coordinates. This still leaves the freedom of 
choice of coordinates within E, which can be redefined, provided this introduces changes 
of the same order as the perturbations 0]. This gauge freedom on E is very important 
because it allows to cast the evolution problem (and corresponding initial data) in different 
physically equivalent forms, that simplify either the mathematical treatment or the physical 
interpretation 

However, in considering more general coordinate transformations one is faced with the 
following problem. Suppose that for certain choice of coordinates x^ = (t,r,9,(p), E is the 
hypersurface t = 0, and that the hypersurfaces E t , corresponding to constant values of t 
define a foliation of the space time manifold Ai, as would be appropriate for the evolution 
of initial data on E, using t as the J'time" parameter. We may introduce now a different 
coordinate system, say x M = (t,f,9,<f), such that the constant (t hypersurfaces provide 
a different foliation of A4, but where the transformation between x^ and x M depends on 
the €i in such a way that we get x^ = x M for e,i = 0. It is not immediate that this more 
general case can be treated as that where we consider coordinate changes only within E. 
This is because, by definition, the initial value problem, and corresponding initial data, 
for the coordinates x,^ correspond now to constant t, but the restrictions on the limits of 
the €i = 0, do not necessarily imply that, for finite e^, even if the hypersurfaces E and E, 
corresponding respectively to t — and t = intersect at some points, they cannot become 
widely separated in time as we move away from these points. (See below for more details). 
But this would imply, in principle, that in this case the only way we can implement the 
coordinate transformation is if we have already solved the evolution equations, since a finite 
amount of time is required to move in general from a point in E to a point in E. 

The previous discussion is relevant to the following problem. Suppose we are given a 
family of initial data for the Einstein equations, (3-metrics gij, extrinsic curvatures ify), 
depending smoothly on a parameter P, such that for P = we recover Schwarzschild's 
space time, and for P ^ the initial data is asymptotically flat, with ADM momentum P. 
Quite generally, for small P, we expect this data to evolve into a "boosted" Schwarzschild 
black hole, possibly accompanied with the emission of a certain amount of gravitational 
radiation. However, if we consider applying the Regge- Wheeler-Zerilli perturbation theory 
to analyze this evolution, we are faced with the problem that this theory is based on the 
assumption that the evolution leads to a stationary black hole, essentially centered at the 
origin of spatial coordinates, while the "natural" evolution of the given initial data leads 
to a non stationary final state of the black hole. Clearly, this problem could be solved 
by choosing a new foliation, where the final state of the black hole is stationary, but such 
transformation, for any P^O involves a "Lorentz boost" , with arbitrarily large separations 
of the hypersurfaces E and E, as can be seen by considering the simpler similar case in 
Minkowski spacetime. 

The evolution of conformally flat "boosted" single black hole initial data in black hole 
perturbation theory, was analyzed in a recent paper by Khanna, Gleiser and Pullin [?}, 
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(referred as I in what follows), but there the previous problem was given only a heuristic 
treatment. In more detail, the analysis performed in I makes use of second order perturbation 
theory, not as regards evolution, but rather to carry out a second order gauge transformation 
that eliminates the first order terms, leaving only the second order contributions, which then 
satisfy linearized Einstein equations. Although the results obtained regarding radiate wave 
forms are qualitatively in agreement with other results obtained in perturbation theory, i.e., 
they show, for instance, the expected "quasi-normal ringing", an intriguing feature that 
distinguishes this case is that the Zerilli function ip does not vanish for large radial distance, 
approaching instead a constant non vanishing value. Since the radiated energy depends on 
the time derivative of ip, the presence of this constant does not in itself mean that there is 
a divergence, but this behavior is in clear contrast with that previously observed in other 
applications of perturbation theory, and, therefore, it justifies a more detailed analysis and 
interpretation. 

The other point that also needs consideration in detail, for the reasons mentioned above, is 
the type of gauge transformation performed in I on the initial data. Its effect was equivalent 
to a coordinate transformation where one moves from a slice where the black hole has non 
vanishing linear momentum, to a frame where it is "at rest", that is, the transformation 
is essentially a "boost". But this introduces, at least in principle, a transformation that 
requires the knowledge of the evolution of the initial data from the "boosted" frame to the 
"rest" frame, and therefore, as indicated, it is not simply equivalent to a relabeling of points 
on the initial data surface. 

In this paper we consider again the problem from a more general point of view. We first 
present a justification for the validity of the "passage to the center of mass system" used in 
jjj, and then show that the somewhat unexpected asymptotic behavior of the Zerilli function 
found there can be interpreted as a gravitational memory effect, that might be present in 
some form in any problem where one has "single boosted black hole" type of initial data. 



II. A DIGRESSION ON COORDINATE AND GAUGE TRANSFORMATION 

The development of higher order perturbation theory given in (see also [H) is based 
on the existence of a family of solutions of Einstein's equations, depending on the parameter 
e, which includes the Schwarzschild metric for e = 0, and on the possibility of performing 
general coordinate transformations, which may also be classified in one parameter families, 
with the same parameter, e, as the family of metrics. It is then assumed that both the 
metric coefficients, and coordinate transformation functions may be expanded in powers of 
e, around e = 0, which naturally leads to a classification in "orders", in accordance with the 
corresponding power of e in the expansions. With these assumptions one obtains an infinite 
set of relations between the metric coefficients corresponding to the same geometrical metric, 
expanded in powers of e, but written in different coordinate systems, and the expansion 
coefficients of the coordinate transformation functions, each member of the set corresponding 
to a given order in e. We generally call "nth order gauge transformation" the relations 
obtained equating coefficients of nth order in e. 

For the purpose of applications of perturbation theory, it is only practical to consider the 
lowest orders. In particular, we may consider coordinate transformations that contain only 
linear terms in e. These naturally, generate first order gauge transformations, which are 
linear in the coordinate ("gauge") transformation functions, but they also generate higher 
order gauge transformations, through terms that are quadratic, cubic, etc, in the gauge 
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functions. Similarly, we may consider coordinate transformations that are quadratic in e, to 
start with. These generate second, fourth, etc, order gauge transformations, but do not affect 
the first order terms. Thus we may consider, as in |5|], a sequence of gauge transformations, 
where the order of the gauge transformation is raised as we move along the sequence. 

In all these considerations, we are assuming that the metric is known in some four dimen- 
sional region of the spacetime manifold, so that the coordinate transformations are quite 
general. We remark, however, that an important set of applications to black hole physics 
is based on the perturbative solution of an initial value problem. This requires the intro- 
duction of some foliation of spacetime that singles out a one parameter family of spacelike 
hypersurfaces, on one of which the initial data is given. The simplest way of specifying 
the hypersurfaces is by introducing a "time" coordinate t, such that t = constant, on each 
hypersurface, and the initial data is given for t = 0. Under a general gauge transformation 
of the kind described above, we might introduce a new "time" coordinate t', such that it 
also provides a foliation of spacetime, and the initial value problem could, in principle, be 
solved starting with initial values on the hypersurface if = 0. 

We notice, however, that the initial value problem, from its formulation, implies in prac- 
tical applications that the metric (and its first t derivative) is known only for t = 0, on 
some "given" hypersurface. If this is all the knowledge of the metric that we have at the 
beginning, and, if the hypersurfaces t — 0, and t' = 0, do not coincide, to obtain the corre- 
sponding initial value for t' = 0, we need in principle to solve the evolution equations, since 
points on t' = 0, may be arbitrarily far to the future, or past, of points on t = 0. Thus, in 
practice, where an initial value problem is concerned, one cannot apply the full set of gauge 
transformations, but a restriction must be imposed so that either, there is no change in the 
initial data hypersurface, or the change is only of the order of e considered. 



A. A toy model 



We may illustrate these points with a "toy model" . Consider a field theory in one-plus-one 
dimensions, where the field satisfies the Sine-Gordon equation, 

d 2 (j){x,t) d 2 <p{x,t) . 

—W ^ + sm0(x,t)=O (1) 

Any solution (f)(x,t) of (JH) is completely determined by the "initial data", 

<f>{x,t = 0) , {d<P{x,t)/dt)\ t=0 . 
Equation (JIJ is invariant under transformations ( "boosts" ) of the form, 

x = cosh(i])x' — sinh(?7)t / 

t = cosh(r/)t' — sinh(^)x' (2) 

if 4>(x, t) transforms as a scalar, i.e., 4>'(x', t') = cj)(x(x', t'),t(x', t')). Moreover, corresponding 
to any solution (j)(x,t) of (JTJ), we may define the quantities, 

f - / \(<P,t) 2 + (<P, x ) 2 + sm 2 <P}dx 



((j),t)((j>,x)dx (3) 
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where the integrals are computed for constant t, but E, and P, are actually independent of 
t, and transform as the (t, x) components of a two- vector, i.e., under the transformation (J2J), 
we have, 

P = cosh (rj)P' — sinh(i])E' 

E = cosh(r])E' - sinh(r/)P' (4) 

where E', P' are related to as E, P are related to <fi. We may think of E and P, as the 
"observable" quantities to be computed from 0. The crucial point here is that E and P may 
be computed solely in terms of the initial data for t = 0, and similarly E' and P' may be 
computed in terms of the initial data for t' = 0. 

To make contact with the discussion in this paper, we may view the transformations (0) 
as defining a one parameter set of solutions of (JIJ, namely, we define 

4>(x, t, rf) = (p (cosh(ri)x — sinh(^)t, cosh(^)t — sinh(77)a;) (5) 

where <f) (x,t) is some solution of (JTJ. Correspondingly, we have a one parameter set of 
2-vectors (E(rj), P(j])), obtained from (p(x,t,rj) through 

Suppose now that we want to consider the evolution equation (JTJ) for the different 4>(x, t, rj) 
as an initial value problem. If the solution <p (x,t) is given, (essentially for all t), then we 
have, 

(f>(x,0,r)) = o (cosh(?7)x, — smh(rj)x) 
<p(x,0,T)), t = -(f) o (x,0), x smh(r)) +<f) o (x,0), t cosh(r)) (6) 

and, somewhat trivially, we can solve the initial value problem for the data 
(4>(x, 0, rj), <j)(x, 0, 7)), t ), to recover <fi(x,t,r)), and (E(r)), P(r))). 

On the other hand, if all that is known for <f) (x,t) is the initial data for t = 0, then, 
we cannot write the right hand side of (0), and we need to solve the initial value problem 
for <ft (x,t) before we can proceed. However, if we are only interested in an "infinitesimal" 
boost, i.e., the limit rj — > 0, we may attempt a computation of the right hand sides of (0) by 
expanding (j) Q in a power series in 7]. Namely, since, 

x = (I+7772 + ...)x' - (77 + ...)*' 

t = (l+r ] 2 /2 + ...)t'-( V + ...)x' (7) 

we have 

(f>(x,0,r]) = o ((l + r/ 2 /2 + ...)x,-(r] + ...)x) 

~ o (x,O) - o (x, 0), t 1] X + ... (8) 

and a similar expression for 0(x, 0, rj), t- But, due to the presence of the factors x on the right 
hand, this means that by restricting to the lowest powers of 77, we may obtain an expression 
for the initial data for (p(x, t, rj) that differs drastically from the exact form. We notice, for 
instance, that while the exact initial data might be square integrable, this might not hold for 
the "perturbative" expression in the right hand side of (jHI)- The problem here may be traced 
to the non commutativity of the limits \x'\ 00, and 77 — > 0. In the "boost" interpretation, 
for any 77 7^ 0, the "hypersurfaces" t — 0, t' — 0, become arbitrarily separated at large \x'\ 
(or \x\), and the evolution equations must be satisfied to move from one to the other. 
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This does not mean that the "boost" transformation cannot be used in a "perturbative" 
sense in any Lorentz invariant model. Consider, instead of the general form for <p (x,t) the 
"solitons" , 

x{x, t, (3) = 4 arctan(exp(cosh(/?)x — sinh(/3)£)) (9) 

where (3 is a constant. 

Let us take (f> Q (x,t) = x( x ,t,(3 = 0). This solution is static, i.e., independent of t, and 
the solutions for f3 ^ are obtained applying a boost with r] = (3 to (f> (x,t). In particular, 
the initial data for <f)(x,t,rj) will be related to <p (x,t) by 

0(x,O,ry) = o (cosh(/3) x, 0) 

= 4 arctan(exp(cosh(r/)a;)) 

(f)(x,0,r]), t = -0 o (cosh(/3) x,0), x sinh(/3)/ cosh(/3) 

= — 2 sinh(^)/ cosh(cosh(r]) x) (10) 

and we notice that the "initial data" for <p Q (x, t) is enough to compute that for <p(x, t, rj), for 
all rj. 

If we consider again an "infinitesimal" boost (77 — » 0) of the static soliton, we find 

Atv 

<b' ~ 4arctan(Aexp(cosh(x)) -4— — ; + 0(v 2 ) (11) 

v v v " A 2 exp(x) +exp(-x) v ' K ' 

and we notice that in this case the "perturbation" is uniformly bounded in x, and we may, 
for instance, use this expression for <p' to compute, say, E or P, to the corresponding order 
in rj. This behavior is quite different from that found for a general non static solution, and 
originates in the fact that the "unperturbed" solution is static, (independent of t). If we 
repeat the arguments for the failure of the expansion in the general case, we notice that for 
the static solution, the initial data is the same as the solution for all t, and therefore the 
initial data on a "boosted" slice is simply obtained by "boosting", i.e. changing coordinates 
in the initial data on the static slice. 

Going back to the black hole perturbation problem, we notice that the Schwarzschild 
metric, in appropriate coordinates, is manifestly static, and, precisely for the same reason 
as above, a coordinate transformation equivalent to a "boost", when applied to this metric, 
although changing the initial data slice, is equivalent to a gauge transformation that does 
not change the slice, because the evolution equations are automatically satisfied. This will 
be considered in detail in the following Sections. 

III. REVIEW OF £ = 1 EVEN PARITY LINEAR PERTURBATIONS 

It will be useful to review some properties of i = 1 even parity linear perturbations of 
a Schwarzschild black hole 0, 0, 0|. Restricting to axisymmetry, these may be written in 
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general in the form, 



» (1 > 

Hrr 

c 

(1) 



(1 -2M/r)iJ ( 5 1) (t,r)cos(^) 

H^\t,r) cos(#) 

1/(1 -2M/r)H^\t,r) cos(fl) 

-^(^r) sin(0) 

-4 1} (t,r) sin(0) 

r 2 K (1) (t,r)cos(0) 

r sin (0)X (1) (t,r)cos(fl) 



(12) 



One can show that the general solution of the (linearized) vacuum Einstein equations 
satisfied by these perturbations can be written in the form, 



H?\t,r 
Hi 1] (t,r 

K w (t,r 



r(r - 2M) 1 v ' ; <9t V ' ; 
r-2M) dt 1 ( ' j r <9r ^ 



r 2 |M«(t,r)-^ 



2M) 



77zW Ml(1)(t ' r)+r2( ¥ M ^ (t ' r)) 



M (1) (*,r) 



M 1 (1) (t 1 r)-2M 2 w (f 1 r) 



(13) 



where and M^ 1 "* are arbitrary functions of (t,r). 

Conversely, given any solution of the (linearized) vacuum Einstein equations, the func- 
tions M- 1 ^ are given by, 



M« = 



Mj 



(i) 



(r - 2M) 
r-2M 



6M 



K^, r -rH^ ] + 2h[ 



,(!) 



(i) 



- V 1} + 



r -2M 



^), r -riJ 2 ( 1 ) + 2M 1) ) 



2"" ' 6rM V " " ■ •-■■!/ ( 14 ) 

An important property of the I = 1 even parity perturbations is their relation to the 
ADM linear momentum associated to an asymptotically flat initial data set (3-metric h^ v 
and extrinsic curvature K^ u ) given on an asymptotically flat hypersurface E. To obtain ex- 
pressions more directly related to black hole perturbation theory, we consider asymptotically 
Euclidean (cartesian) coordinates {x 1 }, and a related set of spherical coordinates {r,9,ip}, 
which we indicate by {x^}. One can show that the Euclidean components of p M are given 
by 



Pi 



8tt 



lim y 



dx v 
dx l 



{K^n* - K%n u ) dA 



(15) 
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where refer to Euclidean (cartesian) coordinates and fi,u... to spherical coordinates. 

The integral is taken over a sphere of radius r, with n M the unit normal to the sphere and 
dA = r 2 sin 9 d9 dip , the area element on the same sphere. 
We recall that if we write the 4-metric in the form 



ds 2 



-N 2 dt 2 + h %j [dx l + N i dt){dx i + N j dt) 



(16) 



where x l are coordinates on the hypersurface, and N and N % are respectively the lapse 
function and shift vector, we have 



dh i3 
dt 



Ni\j + K 



j\> 



2NK, 



(17) 



where "|" indicates covariant derivative with respect to 3-metric gij, on the hypersurface S. 
Then, it is easy to show that if the metric is given in the form (|T2*|) . and ( is a perturbation 
parameter upon which the I = 1 even perturbation depend linearly, equation (fTo*)) takes the 
form, 



r 0" 

hm — 

r— >oo 6 



(i) 



dh\ 



(i) 



dK^ 



dt 



(18) 



We remark, more generally, that only the 
asymptotically to P z . 



dt dr 
1, even parity perturbations contribute 



IV. THE EFFECT OF A BOOST ON THE ZERILLI FUNCTION 

Suppose we have a black hole perturbation problem where the leading perturbations, 
which we take as of order e, (where e is some parameter), are of even parity with L = 2, 
m = 0. Assume further that the metric is written in an asymptotically flat gauge, and 
that first order perturbation theory is appropriate for an analysis of the evolution of the 
perturbations. This implies that both the linear and angular momentum vanish to order e. 
We may construct the corresponding (£ = 2, m = 0, even parity) Zerilli-Moncrief function 
which may be written as, 



2M) / ( 2) _ rR{2) ^ +3 rGf(2) ^ J_ h {2)^ + K (2) 



(19) 



(2r + 3M) 

and the vacuum Einstein equations imply that ip satisfies the Zerilli equation, 

d 2 ib d 2 ib , , 

w-g$- v * = ° < 20 > 

where r* = r + 2Mln(r — 2M), and, for even parity £ = 2, the "potential" V is given by, 

/ 2M\ (4r 3 + 4Mr 2 + 6rM 2 + 3M 3 ) 
V[r) ~ 6 V~~) r3(2r + 3M) 2 (21) 

Suppose further that on some constant t slice we have ip — ► 0, and (dijj/dt) — > as 
r — > oo, with ijj and (dijj/dt) some smooth bounded functions for r > 2M. Then using (|2*U|) 
we generally find that, after an eventual "quasi- normal ringing" type of waveform, ip — > 0, 
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for all fixed r >> 2M, as t — > oo. This in turn implies that, given the simple relationship 
between if) and the gravitational wave amplitude, an interferometric detector, initially in a 
certain state, will return to that state after the passage of the quasinormal ringing wave, with 
no "memory" effect. Examples of this type can be found, e.g., in We also notice that 
generally, if the initial data contains only even L terms, no odd L values will be generated 
in the evolution. 

Let us assume now that e > 0. We may set e = P 2 , and consider P as a perturbation 
parameter. Then, in the spirit of |3], we assume that a full solution of Einstein's equations 
exists, given by a metric # M i,(r, 9, 0, t), for which if) is obtained by an appropriate expansion 
in P. Moreover, we may assume that this metric is given in a coordinate system where it 
is explicitly in asymptotically flat form, and that for P = we recover the Schwarzschild 
form for the metric. We may now consider applying on a coordinate transformation of 
the form, 

t = t' + PM (t', r') cos(0') 
r = r' + PM 1 (t',r')cos(9') 
= 0' - PM 2 (t',r') sin(e') 

= 0' (22) 

Then, (recall the previous discussion) if we expand the transformed metric in powers 
of P near P = 0, we find that the zeroth order terms take again the static Schwarzschild 
metric form in (r' , 9' , <f> , if) coordinates, while to order P we find £ — 1, even parity terms 
of the form (|12j) . with the coefficients given in terms the Aii by expression of the form 
(I13j) . where the functions Mj are replaced by Aii. This is quite general, but, to simplify 
the physical interpretation, we may now restrict the functions Aii in such a way that the 
resulting i = 1, even parity terms satisfy the conditions of asymptotic flatness. By this we 
mean that asymptotically for large r' we should have (t , r') , H^\t',r') and K^(t',r') 
at most of order 1/r' 2 , and H^\t',r'), h^\t',r') and h\ (t',r r ) at most of order 1/r'. One 
can then check that this is possible only if the functions Aii are asymptotically of the form, 

M$ = a r' + ai + a 2 /r' + 0(l/r' 2 ) 
M ( P = a t' + a 3 + 0(1/0 

M { 2 1} = (a t' + a 3 )/r' + ((ai - 2Ma )t' + a 4 )/r' 2 + 0(l/r' 3 ) (23) 

where the are constants. Moreover, if we replace the resulting expressions for flf^f.O* 
h { o\t',r'), h^\t\r') and K^\t' , r') in (HHJ) we find, 

P z = lim P (Ma + 0(1/0] ( 24 ) 

r— >oo 

Therefore, if we choose, 

a = 1/M, (25) 

the I = 1 even parity part of the metric describes an ADM momentum equal to P, and the 
transformation (|22j) carries the metric to a "boosted frame", with momentum P. 

The preceding discussion is nevertheless incomplete. The reason is that actually equation 
()18|) and its physical interpretation hold only on an asymptotically flat frame, and one can 
check that the transformation (J2~2")) . even with the restrictions ()23|) . (|23|) introduces terms 
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of order P 2 in the I = and 1 = 2 even parity part of the metric that are not compatible 
with (explicit) asymptotic flatness. This, however, does not modify the interpretation of 
P as the ADM momentum. The reason is that we can restore explicit asymptotic flatness 
introducing a new coordinate (gauge) transformation, of order P 2 , that involves only the 
I = and I = 2 even parity terms. This has no effect on the even parity, i = 1 terms and, 
therefore, leaves the right hand side of (fT%|) unchanged. 

Summarizing, we see that to change from the initial asymptotically flat slicing with 
vanishing ADM linear momentum to a new asymptotically flat slicing where this takes the 
value P, we need to perform at least the equivalent of a gauge transformation of order P, 
followed by one of order P 2 that has no effect on the ADM momentum. 

But now we are in position to discuss the relationship between the Zerilli-Moncrief func- 
tion if) given by (JT5Jl in the zero momentum frame, and the corresponding quantity ips 
constructed using equation (fTTIJ) . with the functions on the right hand side as given in the 
momentum P frame. These two functions are not equal but, considering now the effect to 
order P 2 of the transformation (|22|) . (i.e., as a second order gauge transformation) one can 
verify that ips — i> can be written as a quadratic homogeneous expression involving only 
the M.i and their t' and r' derivatives. The explicit expressions are rather lengthy, but the 
important result is that with the restrictions (|23|). (|23j) we find that for large r', 

2P 2 (17- 16 ai -a\)P 2 . /2 , 

1">-1> = - m + { ~ Q ;, + 0{iy 2 ) (26) 

Therefore, for large r and any finite t, ips — ip approaches the constant value —2P 2 /3M 
irrespective of the details of the 1 = 2 data. The crucial point that makes this result non 
trivial is that ipB — if> is not changed by the subsequent transformation of order P 2 that 
restores asymptotic flatness, because the expressions for both if) and if)g are gauge invariant 
under those transformations. 

An immediate consequence of these results for initial data sets {gBij, Kb^} with non 
vanishing linear ADM momentum P, is that if the function ips vanishes for large r, then 
the corresponding Zerilli function in the "center of momentum" frame (i.e. the frame with 
vanishing ADM momentum) will approach the constant value 2P 2 /3M for large r. This is 
precisely the behavior observed in I, for the particular case of a conformally flat initial data. 
In the next Section we explore further the properties of this type of perturbations. 



V. A GRAVITATIONAL MEMORY EFFECT 

Consider again the Zerilli equation (J20|) . The properties of the solution obtained by 
evolution of initial data of compact support have been extensively studied following the 
original work of Price and Wald In the case of interest in the present analysis, 

however, the initial data is only assumed to be smooth and uniformly bounded, with dif>/dt 
vanishing for r* — > ±oo, while if) may approach constant non vanishing values as r* — > ±oo, 
and, therefore, the results obtained in the case of compact support are not immediately 
applicable. We may resort, nevertheless, to plausible, albeit non rigorous, arguments to 
predict the evolution of this type of data. As we shall see, the results are in agreement with 
what we obtain by numerical methods. 

First we notice that the "energy" integral 

/+oo 
[(dif)(t, n)/dt) 2 + (dif)(t, r*)/dr*) 2 + if)(t, n) 2 V{n)} air* (27) 
-oo 
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is finite for the initial data (at t = 0) that we are considering, provided only that d t ip — > 
sufficiently fast for r* — > ±oo. 

Moreover, since at large |r*| the Zerilli equation approaches the free wave equation form, 
the limits of ip for r* — > ±oo are not changed by evolution through a finite time. Therefore, 
S(t) should be constant in time, because, on account of (|2()jl . we have, 

dt W r^+oo\dt J \drj r^-oo\dt J \drj v ; 

and therefore d£(t)/dt = 0, for finite time. 

Assume now that M is the black hole mass, and K\ and K 2 are large positive numbers. 
We may consider now initial data T> Q that coincides with our data D inside an interval 
X = (—KiM, K 2 M) of r*, but is of compact support outside of X. We expect that the Zerilli 
function ip resulting from the evolution of T>, in the domain of dependence of X will display 
a "standard" behavior, namely, at sufficiently large t, we expect ipo to display a quasi-normal 
ringing wave form, plus a "tail", and essentially vanish for fixed r#, after the passage of the 
quasi- normal ringing signal. But, since T>o and T> coincide in X, this should also be the 
behavior of the function ip resulting from the evolution of T> in the domain of dependence 
of X. 

We may also get an idea of the behavior of ip at large r, (this is somewhat simpler than 
for r*, because it avoids irrelevant logarithms) if we assume that at t = we have d t ip = 0, 
and tp admits an asymptotic expansion of the form, (as is true, for instance, for the data of 

0), 

^{r, t = 0) ~ C + d/r + C 2 /(r) 2 + ... (29) 

In this case, an asymptotic expansion ip(r,t), satisfying this initial data and (|2*0|) is of the 
form, 

t) ~ C + Ci/r + C 2 /r 2 - 3C t 2 /r 2 + 0(t 2 /r 3 ) (30) 

Then, for large r we may assume that the 0(t 2 /r 3 ) terms are negligible for t < r, and even 
for t ~ r. If this is correct, for sufficiently large t, irrespective of other details, psi should 
change sign around r ~ and smoothly approach the constant value Cq as r — » +00. 
Moreover, we find that d r tp ~ 2Cq/t, and d t ip — 2v / 3Co/r at point where ip changes sign, 
and therefore both approach zero as t increases. Notice that the point where ip = 0, "moves 
faster than light" . This, however, is only a consequence of the form of the initial data, and 
does not imply any causality violation. 

Putting these results together, we expect that after some time to of evolution of our 
initial data, we should see an essentially vanishing wave form up to r ~ to, where we should 
see some quasinormal ringing type wave form, followed by a region where ip slowly changes 
sign and finally approaches the constant asymptotic value of the initial data. This is indeed 
what happens if we numerically integrate the Zerilli equation. In Figure ^ we indicate the 
results of evolving initial data of the form ip{ r *, t = 0) = 1, d t ip(r*, t = 0) = 0, form M = 1, 
through several times. Curve a) corresponds to t — 520M, curve b) to t = 1041M, and c) 
to t = 2082M. A detail of the quasinormal ringing part, for t = 510M, is included in Figure 
2, a), while Figure 2, b) displays the time derivative of ip. It can be seen that these results 
display all the expected features. In particular, we find a form of "gravitational memory 
effect", since at large distances from the black hole, the gravitational wave amplitude ijj 
starts with a certain non vanishing constant value and goes to a different constant value (in 
this case zero), after some quasi-normal ringing signal is observed. We discuss this effect in 
more detail in the next Section. 
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FIG. 1: The evolution of initial data of the form {4>(t = 0,r) = l,ipt(t = 0,r) = 0}, at times 
t = 520M, (plot a)), t = 1041M, (plot b)) and t = 2082M, (plot c)). The quasinormal ringing, 
noticeable as a structure near r = 500M, r = 1000M, and r = 2000M, respectively in a), b) and 
c) is shown enlarged in Figure 2, a), for t = 520M. Notice that the point where ip first changes 
sign moves to the right at about twice the speed of the quasinormal ringing, as indicated in the 
text 



VI. FINAL COMMENTS AND CONCLUSIONS 

The results obtained in the previous Section indicate the existence of a certain grav- 
itational memory effect associated to a particular type of initial data, for instance, that 



obtained for single boosted black holes using the Bowen - York Ansatz . We may picture 
what this implies by considering its effect on an interferometric type detector of gravita- 
tional waves, placed at a large distance from the black hole, that is turned on at t = 0. 
There one finds a very slow drift of the equilibrium position, followed by the quasi-normal 
ringing signal, with the interferometer ending up in an equilibrium position that is shifted 
with respect to the initial one, in an effect that holds certain similarity to the well known 
"gravitational memory" effect |l2l . fbi. fl3 | . In fact, this is closer to the type of effect origi- 
nally envisioned by Grischuk, et.al. [l5|, while its relation to the Christodoulou |12| type of 
gravitational memory is unclear to us, although we remark that the nonlinear nature of the 
Einstein equations is involved, since we had to consider second order gauge transformations 
to arrive at our final result. 
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FIG. 2: The region where the wave form is dominated by quasinormal ringing, at t 
indicated for ip in a), and for dip/dt in b) 



520M, is 



It is well known that the initial data constructed in accordance with the conformal flat- 
ness prescription has the physical drawback that it generally implies incomi ng g ravitational 
radiation in the past. For binary black hole collisions, such as in the Misner [16[ initial data 
or, more generally, for boosted black hole data as considered in this seems to introduce 
no particular undesirable or unexpected features, at least in the close limit approximation 
in center of momentum frame. For a perturbed single boosted black hole, however, we find 
that the conformal flatness of the initial data introduces a new type of gravitational memory 
effect that we would expect to be absent in the absence of incoming gravitational radiation. 
Going back to the original problem of the "appropriate" asymptotic behavior of the initial 
data for a perturbed boosted single black hole, we may conclude that the "natural" choice 
that would be physically expected in the absence of substantial incoming radiation, should 
lead in the boosted frame where the momentum is P, to a Zerilli function that approaches 
the constant value ip — > (2/3)P 2 /M for large r , in agreement with the discussion carried 
out in this Paper. 

Finally, our results can also be considered in relation with analyses such those of Ken- 
nefick [l8|, applied, for instance, to a system such as a black hole surrounded by a spherically 
symmetric stationary mass distribution, that eventually collapses in an asymmetrical man- 
ner, inducing a recoil of the final black hole, together with the emission of gravitational 
radiation. In accordance with our analysis, far away from the source, we have initially a 
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stationary Schwarzschild spacetime, where ip vanishes, but, after the passage of the gravita- 
tional radiation associated to the collapse, the spacetime should correspond to a "boosted" 
Schwarzschild, where ^ is a non vanishing constant, exactly as envisioned in 0|. Notice 
that the final value of ijj is related to the recoil momentum as derived in Section XX, and 
this in turn is related to the anisotropy of the emmited radiation, so that we seem to have 
a further interpretation for the Cristodoulou type of gravitational memory. 
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APPENDIX: A BOOSTED SCHWARZSCHILD METRIC 

Consider the Schwarzschild metric, written in the form 

ds 2 = -(1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 [d6 2 + sin 2 ^0 2 ] (A.l) 

We consider (JA.1J) only for r > 2M, and all t. In that region, we may consider, instead of 
t, r, 9, 0, new coordinates t', r', 9 1 , <fi', related to the old ones by 

= 0' 

9 = 9'- (vt'/r') sin(0') 

r = r' — vt' cos(6*') 

t = t' - vr' cos(9') (A.2) 

This coordinate transformation is motivated by considering a Lorentz boost along the pos- 
itive z-axis, with velocity v, of an auxiliary coordinate system t, x, y, z, related to t, r, 9, <f) 
as flat cartesian coordinates to the corresponding Schwarzschild spherical polar coordinates, 
and keeping only terms of order v, but this is not central to our discussion. The transfor- 
mation (jA.2|) can be used to write the Schwarzschild metric in t',r',9',(j)', coordinates. The 
resulting expression for the line element is rather long and we shall not display it here. The 
important point is that it is diffeomorphic to (|A.1|) . and the metric components are smooth 
functions of v near v = 0, with (jA.l|) recovered for v = 0. If we expand the metric coefficients 
in powers of v, we find that to order v the metric is asymptotically flat in the sense used 
in Section IV. We may use (|18|) to compute P. The result is P = vM, which corresponds 
to the relativistic momentum of an object of (rest) mass M, computed to order v. If we 
consider now the expansion to order v 2 , we find that the metric is not asymptotically flat 
to that order, again in agreement with Section IV. We may now consider a new coordinate 
transformation, of order v 2 , to restores the asymptotic flatness, but, since the Zerilli function 
is invariant under this transformation, we may use the results obtained already in t', r', 9 ; , </>' 
coordinates to compute it. The result is, 

A Av 2 M(r' -2M)(3r' + 2M) 2 , , (N , , 

$ B (t' = 0, r') = i — -. — — — -r— ~ — Mv 2 + 0(1 r') A.3 

Y v ' ; 9r' 2r' + 3M 3 y 1 ' y ' 
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which displays the expect asymptotic constant value for ip B ~ 2P 2 /(3M). 
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